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Constraints on cosmological models from strong gravitational lensing statistics are in- 
vestigated. We pay particular attention to the role of the velocity function in the calculation 
of the lensing probability. The velocity function derived from the observed galaxy luminos- 
ity function, which is used in most previous work, is unable to predict the large separation 
lensing events. In this paper, we also use the Press-Schechter theory to construct a velocity 
function theoretically. Model predictions are compared with the observed velocity function 
and the HST snapshot survey. Comparison with the latter observation shows that the pre- 
dictions based on the theoretical velocity function are consistent with the observed large 
separation events in COBE normalized low-density models, especially with a non-vanishing 
cosmological constant. Adopting the COBE normalization, however, we could find no model 
which simultaneously satisfies both the observed velocity function and the HST snapshot sur- 
vey. We systematically investigate various uncertainties in the gravitational lensing statistics 
including finite core radius, the distance formula, magnification bias, and dust obscuration. 
The results are very sensitive to these effects as well as theoretical models for the velocity 
function, implying that current limits on the cosmological parameters should be interpreted 
with caution. Predictions for future surveys are also presented. 



§1. Introduction 

To date, it is well known that the probability of gravitational lensing (hereafter GL1 of high 
redshift objects is sensitive to the cosmological constant A through the geometrical effect 0)'fl '!)'□) 
and that extremely Tl-dominated universes predict too many lensing events which is inconsistent with 
the observed low frequency of lensed quasars in current samples ErEP> IIP. KochanekcP, for example, 
concludes that Ao := A/(3Hq) < 0.66 at 95% confidence. 

GL limits on Ao should be contrasted to the other cosmological tests which favor low density 
universes with non- vanishing Ao- The latter tests include the two point angular correlation function of 
galaxies B, thegalaxy number counts©, the cluster mass and peculiar velocity functions0'©'©'0 
(See Refs.^'Ef for reviews). Moreover, the recent measurements of the (local] Hubble constant 
H =: lOO/ikms^Mpc -1 give somewhat large values around h ~ 0.t0.8tZP'li3)'l!J)'E$. To reconcile 
the age of the universe ~ -f^o" 1 with that of the oldest globular clusters c3 and of the young galaxies EHP, 
it seems inevitable to introduce non- null A of A > 0.8§). Therefore it is worthwhile reexamining the 
previous constraints on A from GL statistics by taking account of several uncertainties involved in the 
calculation of GL probability. 

In this paper, we pay particular attention to the role of the velocity function (hereafter VF) of 
lensing objects. Previous work on GL statistics usually uses VF determined from the observed galaxy 
luminosity function. Moreover, it is often assumed that the comoving number density of lensing 
objects is constant. However, the luminosity function counts only luminous objects while invisible 
dark objects might as well be responsible for GL. Bearing this in mind, we also use the VF derived 
from the Press-Schechter theoryES (hereafter PS) with some specific cosmological models such as the 
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cold dark matter (CDM) and the primeval isocurvature baryon (PIB) models Bj). The theoretical 
velocity function, by its construction, properly counts all gravitationally bound dark halos as well as 
luminous objects, and automatically takes account of the hierarchical evolution of objects for various 
cosmological parameters^. Narayan & White EIP first applied the PS theory to GL statistics, but only 
in the standard CDM cosmology with no comparison with observation. KochanekE^ extended their 
calculation for various cosmological parameters, and argued the constraints on the bias parameter b 
within the standard CDM model. However, we will see that the PS theory in its original form is not 
a perfect one in deriving VF, and adopt theoretical attempts towards more realistic YF. We compare 
these theoretical VF and their predictions on GL statistics with the observed VFc3) and the HST 
snapshot survey® for COBE normalized CDM and PIB models. 

We also make quantitative estimates of following uncertainties in GL statistics: the effect of 
inhomogeneity of the universe in the distance formula, the finite core size of lensing objects, the 
obscuration of images by dusts in the lensing galaxies, andthe di ff erent definitions of the magnification 
bias. Although each topic is already discussed separately^' IP' O'^'EiJ'ES'LS', we take account of all 
the above effects consistently and systematically for the first time. 

This paper is organized as follows. In §|2], we briefly summarize the basic formulae of GL statistics. 
In §|3] we present the VF's which are applied to GL statistics in the later sections. In §Q, we compare 
the theoretical VF and its GL prediction with observation. In §|5|, we survey various uncertainties 
described above and see how the results in §|4] are affected by them. In §||, we give predictions for 
future GL surveys. Finally in §|7], we summarize the main results of this paper. Our main results are 
discussed after §4 and are summarized in Table II. Those who are familiar with GL and the PS theory 
can skip §2 and §3 which are written only for completeness and terminology. 

Throughout the paper, we use the units c = G = 1. The cosmological parameters are defined as 

8np x A ^ 8irp A 
H:=-= H [f2 a- 3 + (1 - fl - X )a- 2 + A ] 1/2 , (1-2) 



Q '- 3fl*' X '-3IP> ^'-JM' Ao: "3M' 



a 

where p = po/a? = /?o(l + z) 3 is the mean density of the universe at redshift z and the subscripts 
mean the present epoch z = 0. 

This paper is largely based on Ref.E^ (unpublished) which contains further detailed discussion. 

§2. Brief summary of GL statistics 

2.1. Lens model and GL cross section 

We assume that the lensing objects are spherical and isothermal in describing GL statistics. 
The spherical lens model gives sufficiently accurate description of the properties of lensed images 
such as the image separation and the magnification, as compared with more realistic elliptical lens 
models The isothermal profile is characterized by the internal one-dimensional velocity 

dispersion vc$. Thereby we start from the lens potentialclF (twice the Newtonian potential integrated 
along line-of-sight): 

m = ^v 2 (e+e c ) 1/2 , (2-i) 

where £ istfie impact parameter in the lens plane and the core radius £ c eliminates the central 
singularity.^ The lens model with £s-,= is the "singular isothermal sphere" (SIS). Due to the 
spherical symmetry, the lens equation eIP becomes a scalar equation: 

V=^Z-D LS a(0, (2-2) 



*^ The parameter £ c defined here approximately corresponds to twice the core radius of Ref.l 
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where 77 is the source position in the source plane, -Dol 5 -Dos and Dls are the angular diameter 
distances between the observer, lens and source, and a(£) := (d/d£)ip(l 
thin lens approximation, which is implicit in Eq.( |2-2| ), can be justified c 
GL, because highly virialized and compact galaxies or clusters are considered here as the lensing 
objects. Defining 

,.2-Dol-Dls „ ._ „„„. 2; 



is the deflection angle. The 
•Ef'cS in the case of strong 



Attv 



x :- 



Dos 

y-=v/v* 



4irv D LS , 



one can rewrite Eq.( |2-2j ) as 



x{x 2 + xl)- 1 ' 2 . 



(2-3) 
(2-4) 

(2-5) 



Given a source position y, the image position x is obtained by solving Eq.(2-5). If x c < 1 and 
y < y T := (1 — x^ 3 ) 3 ^ 2 0(l — x c ), triple images form [0(x) is the step function]. The brightness of an 
image at x is magnified by a factored H P (x) := \(y/x)(dy/dx)\~ 1 . When x c / 0, we solve Eq.(^) 
numerically for the image positions Xi(y) {i = 1,2,3; x% < 22 < £3) and calculate Ac(y) := X3 — x\ 
(separation of the outer two images), fj,(y) := YUthfei) (total magnification of all the images), and 
r(y) := exp | ln[/x p (x3)//j p (xi)]| (brightness ratio [> 1] of the outer two images). As noted in Ref. 
it is safe to approximate Ac(y) ~ Zte(0) = 2x% := 2(1 
radian is 



8nv 2 S^(l 



Xc) 1//2 (9(1 — x c ). The separation angle 9 in 
xlf'\ (2-6) 



£>OL -Dos 

The GL cross section© o~(Q) (defined as the area of such a region in the source plane that a source 
is lensed with certain properties Q if it resides within the region) is generally written as the integral 
in the y-pktne: 

poo 

&(Q) = 2tt V 2 S(y,Q)ydy. (2-7) 
J 

For example, if Q is "formation of triple images", then S(y,Q) = 0{y T — y) and a = 7r{r]*y r ) 2 ; if Q is 
"image separation is larger than Zkc*" and "magnification is larger than /j*" and "brightness ratio is 
smaller than r#", then 

a(Av*, /i*, r*) = rjl 0(2x t - Ac*)a(^*, r*) , (2-8) 



where 



tr(A**, r*) := 2vr / y - /j*) (r* - r(y)) 



(2-9) 



2.2. Distance formula 

The ^giilar diameter distances (Z?ol 5 -Dos, -Dls) are given by the solutions of the Dyer-Roeder 
equation: 



i6|),^j). 



d 2 3 + q(z) d 3 afl(*) 
dz 2 ' + TT7^ + 2 (1 + z) 2 



D = 0. 



(2-10) 



where q := 4,0 — A is the deceleration parameter, and the "smoothness parameter" a (0 < a < 1) 
measures the degree of inhomogeneity of the universe (for derivation, see Ref.0)). The last term of 
Eq. (|2T0|) , the Ricci focusing term, represents the matter density inside the beam of photons. Denoting 
Da(z*, z) as the distance from z t to z (> z»), the initial conditions of the differential equation ( |2-10| ) 
are 



Da^z*^ z*) — 0, 



d_ 

dz 



D& (-2-* , Z) 



(1 + z*)H{z* 



(2-11) 



In particular, we write D a (z) := Dq(0, z). The solution with a = 1 is the "filled-beam distance," and 



gives the standard angular diameter distance in an exactly homogeneous FRW space-timec£P, while 
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that with a = is the "empty-beam distance." In the derivation of Ec 
focusing due to the tidal shear (Weyl focusing) is negligible. It is shown t 



it is assumed that the 
* that this assumption 



is plausible in most astronomical situations. In appendix [A|, we present specific solutions of Eq.( 2-10| ). 
2.3. GL probability formula 

The probability that a source at redshift z s is lensed with properties Q is generally written as0) 



P{Q,z s 



l 



dz 



(1 + *) 
H(z] 



-D\{z) 



dxv(Q,X,z,z s )N x (x,z) 



(2-12) 



where the integral variables z and x are > respectively, the redshift and the parameter of the lensing ob- 
jects, with the latter completely characterizing their properties (x is not necessarily a one dimensional 
quantity), and N x (x, z)dx is the function" , i.e., the comoving number density of lensing objects 
with the parameter \ ~ X + dx at z (the reason why the filled-beam distances are used in Eq.[ [2-12 | is 
discussed in §5.1). 

In the above lens model (Eq. | ]2-1| ) , the lensing cross section is characterized by a single parameter 
v (assuming that £ c is empirically related to v as Eq. [|2T6|1 ). This is why VF plays a fundamental 
role in the calculation of GL probability. Substituting Eq.([2-£|) into Eq.(2T2), the probability that 
a source at z s is multiple-imaged with image separation angle larger than 6 and total magnification 
larger than and the brightness ratio smaller than r* is 



P(Q) fi* , r* , z$ 



16tH 



dz 



[1 + z) 
H(z) 



Dl(z)D\(z,z s ) \ dvv 4 N v (v,z)o-(p*,r :¥ ) 



(2-13) 



The limits v\ and vi (> v\) of the u-integral are determined by Eq.( |2-6D (see Eq. |B-1[| ). If x c {v) does 
not increase with v (i.e., ^ c increases less rapidly than oc v 2 ), v\ is given by the sole positive solution 
of Eq.( |2-6| ) for given 9, and ^2 = 00 formally (Strictly speaking, setting v 2 = 00 is inconsistent with 
the isothermal lens model because the isothermal profile does not continue to infinite radii. However, 
one can show that a truncation of isothermal profile at the virial radius makes little difference). If 
x c {v) increases with v, then v\ and V2 are two positive solutions of Eq.( |2-6| ). Probability distribution 
of the image separation angle is 



P e (6,n*,r* 
16vr 2 



dz 



d_ 

~~ ~d0~ 

:i + z? 



H(z) 



P(9,fj,*,r*,z s ) 
Df(z)Dl(z,z s ) 



dv 

1b 



v 4 N v (v,z)a(n*,r*) 



1'2 



Assuming 



dv/d6 in Eq.(2-15) is calculated from Eq.( |2-6| ) as 

v 1 



(2-14) 
(2-15) 

(2-16) 



dv 

de 



x 2 c {v) 



62 



px 2 (v) 



In §5.2, we will see that observations suggest p 
v\ and V2 in the cases p = 0, 1, 2, 3 and 4. 



(2-17) 

3. In appendix |Bl we give analytic expressions for 



2.4. Magnification bias 

Because GL causes a magnification of images and bright QSO's are easy to discover, lensed 
QSO's have relatively high probability of being included in a QSO sample. This selection effect 
(magnification biasEf) is computed as follow. Let ^^(L)dL be the luminosity function of sources, 
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and let ^q(L) := <P ( ^(L')dL' . When_one searches for lensed QSO's of the observed flux brighter 
than S, the GL probability increases asE^: 

1 f°° 

Pf(6,S,r*,z s )-.= —— dfM*P gfH (6,iM*,r*,z s )$q(Lii/n*), (2-18) 
9q{L) Ji 

where Pg^ := — (d/dfi^Pg, /2 := [Da{z s ) j Di(z s )] 2 , and the luminosity 

L := 4vr(l + z s ) 4 Dj(z s )(l + z^ 1 S (2-19) 

must be calculated from S in the same way as in the determination of <^(L). The last factor (l+z) 7_1 
is the JT-correction , which assumes that the energy spectrum of QSO's is of the form E oc v" 1 '. 
We put the factor p in Eq.( |2T8| ) because the 'jimlensed" sources are also magnified (if a < 1) by 
that factor on average from the flux conservations. Substituting Eq.(2T5) into Eq.( [2T8] ), the biased 
probability Pf is expressed as Eq.(2-15) with a([i*,r*) replaced by (cf. Eq, f2-9fl ): 



Similarly, for a QSO of observed flux S, the GL probability is given by Eq.(2-15) with cr(/i*, r*) replaced 

^-it)l> yye{r -- r{y)) <^m (2 ' 21) 

The observed QSO luminosity function is fitted by the two power-law mo del@ 

<P®(L)dL = fa[(L/L*) a + (L I L^fY x {dLj L* ) (2-22) 
with the luminosity evolution 

L*(z s ) = L*(l + z s ) k \ (2-23) 

where a = 3.79, = 1.44, fa = 6.4 x 10~ 6 /i 3 Mpc~ 3 , k L = 3.15, and the absolute B-band magnitude 
corresponding to Lq is Mq = —20.91 + 51og/i. These values are obtained from z s < 2.2 QSO sample, 
assuming (f2o, Ao) = (1, 0) and 7 = 0.5 in Eq.( |2vL9| ) (but the values of a and (3 do not depend on the 
assumed cosmological models). Following Ref. EiZP, we assume that Eq.( |2-23| ) is valid up to z s = 3 and 
that, for z s > 3, L*(z s ) = L^3.2^~^^ a -^ and fa{z s ) = ^.Tr^ 11 ^^- 1 )/^-^ with the same 
values of the parameters. The behavior of this luminosity function at z s > 3 is such that the bright 
end keeps constant while the faint end slidesdown by a factor of 3.2 in every unit redshift. We also 
assume that the QSO sample in Boyle et al.E^P itself does not suffer from the magnification bias, and 
neglect the possible magnification due to microlensing by stars in the lensing galaxies Ea). 

§3. Velocity function 

3.1. Schechter VP: observational viewpoint 

Standard calculation of GL probability in most of the previous work uses VF from the observed 
galaxy luminosity function of Schechter form at z = 0: 

Nl ch {L)dL = ML/L*)- a e- L/L *dL/L* . (3-1) 



*-* we use (y r — j/)' 3 /3 ^ 2 as an integral variable in Eqs.( |2-2l[ ) and ( 2-20[ ) in order to eliminate the divergence of the 
integrand on the radial caustics y = j/ r 
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Autofib survey© concludes that a = 1.09, 0* = 0.026/i 3 Mpc -3 and the absolute magnitude corre- 
sponding to L* is M* = —19.20 + 5 log h in B-band. Combining with the Tully-Fisher or Faber-Jackson 
relation L/L* = (v/v*) 7 , Eq.Q3-l|) yields 



N* ch (v,0)dv 



i=E,S0,S 



v -7j(a-l) 


" / v y*- 


dv 


J exp 










V 



(3-2) 



2!)) 



(hereafter SchVF). Therein constant comoving number density of lensing objects, i.e., N v (v,z]_ 
N v (v,0), is often assumed. Using the B-band Tully-Fisher and Faber-Jackson relations in Ref. 
one finds (7,u>) = (2.9, 126 kms -1 ) for S galaxies, and (7, v*) = (3.3, 175 kms -1 ) for E and SO galaxies 
at the above value of L*. The morphological composition is 0*e + </>*so = 0.44^* and <^>*s = 0.56(/)*^ ■ 
We adopt the above values. 

Several problems in applying SchVF to GL statistics include: (1) SchVF counts only luminous 
objects though they are not the only lenses: invisible dark haloes might as well be responsible for GL. 
(2) It is possible that the assumption of no evolution makes the SchVF Jatrongly affected by the local 
and recent property of the universe around us. In fact, Autofib survey E!P found that the luminosity 
function steepens with redshift up to z = 0.75 at the faint end. Some authorsE^'LlP'E^'Bj) consider 
specific evolution models, but their conclusions does not seem decisive in that the models include some 
arbitrary free parameters. [Inclusion of dark objects (1) and the evolution effect (2) into GL statistics 
would, in principle, lead to tighter limits on A.] (3) The Tully-Fisher and Faber-Jackson relations 
may not be universal in particular at high z. (4) SchVF has an uncertainty in the velocity dispersion 
of early-type galaxies by a factor 1.50), which changes the GL probability by (1.5) 2 . [See Ref.& for 
a detailed study against these arguments.] (5)^lusters of galaxies, which are not counted in SchVF, 
can also contribute to strong lensing of QSO'sB'E^. Because of these uncertainties, we would like to 
introduce another approach below to the computation of VF. 



3.2. Press- Schechter VF: theoretical viewpoint 

Narayan & White© and Kochanekll) applied the PS mass function (see, e.g. 
duction) to compute VF: 



Ref.EP for intro- 



Nff(M, z)dM 



2 po 5 c0 (z) 
ttM a(R) 



diner 



dhxM 



exp 



&(*) 
2a 2 (R) 



dM 



where 



dk k 2 P(k)W 2 (kR) 



(3-3) 



(3-4) 



is the rms of linear density fluctuation today on the comoving scale R = [2M/(J?o#o)] 1/3 [P(fc) and 
W(kR) are the power spectrum at z = and the fc-space window function] and 5 c o(z) is the critical 
density contrast extrapolated linearly to today in order to virialize by z (Eq.[C-30]). Assuming the 
isothermality of the lensing objects again, the one-dimensional velocity dispersion is related to the 
mass through 



M 



2r v 



1/2 



-[Qo^iz^HoR 



(3-5) 



where r v is the virial radius of the object and $vo(z) := (R/r v ) 3 is the overdensity of an object at 
z = which virialized at z (Eq.| C-29|] ). Calculation of 5 c q and i9 v o in A = universe is found in 
Ref. t3 . We give analytic formulae for them in i? + A = 1 universe in appendix [C], using the spherical 
collapse model. We normalize a(R) so that o"(8/i _1 Mpc) = (6 is the bias parameter) and use the 
top-hat window function W(x) = 3(sinx — xcosx)/x 3 . From Eqs,( |3-3|) and ( [3-5D , VF is constructed 
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theoretically asEa^'N) (hereafter PSVF): 



N? s (v,z)dv 



1 8co(z) 



(2vr) 3 /2 R3 a (R) 



din a 



dlnR 



exp 



&(*) 
2a 2 (R) 



dv 
v 



(3-6) 



with the v-R relation in Eq. fl3-5|) , PSVF, by its construction, properly counts all gravitationally 
bound objects (including the clusters of galaxies) as well as luminous objects, and automatically takes 
account of the hierarchical evolution of objects in the universe for various cosmological parameters!^. 
The theoretical foundation of theJRS_mass function is now more secureE3) and it is also supported by 
cosmological iV-body simulationsE^'Ea). 

PSVF strongly increases with Qq and weakly dependent on Ao- The former is because the comoving 
scale R is smaller for larger Qq for fixed v (Eq.[3-5|): since the density fluctuation is large on small 
scales [o~(R) is large], larger f2 universes have more objects at fixed v. The latter is because 5 c q and 
$ v o do not so much depend on Ao at low redshifts. At higher redshifts, however, <5 c o and i? v o are more 
dependent on Ao and so is PSVF. The bias parameter b affects PSVF most strongly: since cr(R) oc 
PSVF behaves like oc 6exp(— b 2 ). The dependence on h and 4?b is very weak because they enter only 
through r (see Eq.f4d|) 



3.3. Formation rate in theoretical VF 

PSVF is derived on the assumption that all objects have just virialized at the moment when one 
evaluates VFL2)'lJ§ (see the definition of "i? v o)- I n the application to GL statistics here, this means 
that the deflection of light and the formation (virialization) of the deflector occurs simultaneously. 
In realitv, however, a virialized object survives for some time until it is merged into more massive 
objects E£rEP, and hence there must be an accumulation of objects in a specific velocity range which 
formed at some earlier epochs. Therefore, we must take account of the formation epoch distribution 
of objects in the universe to avoid the above assumption. 

Kitayama & SutoE3 proposed a practical prescription in obtaining the formation epoch distri- 
bution from the random- walk method ES', and calculated theoretically the comoving number density 
F(M, Zf, z)dMdz{ of those objects which formed at Z{ ~ Zf + dzf with mass M ~ M + dM and survive 
without destructed (absorbed into larger hierarchies) until z. Integrating over Zf, more realistic VF 
may be written as (hereafter KSVF): 



N* s (v, z)dv = dv I dz { F(M, z { , z) 



dM 
dv 



(3-7) 



The explicit form of F is found in their Eq.(23). Compared with PSVF, KSVF is shifted to the 
direction of larger velocity relative to PSVF, especially on small scales (see Figs.l and 2). This is due 
to the accumulation of objects, as mentioned above: because smaller objects (at z = 0) are likely to 
have formed earlier in the bottom up scenario, and because older objects have higher velocities (see 
Eq. |3-5| ), the velocity of small objects becomes larger when we take the formation history into account. 



Eq.(3-7) refines the or igin al PSVF in the sense that the formation epoch is taken into account, though 
not fully satisfactorilyLiiP. 



§4. Comparison with observation 



4.1. Comparison with observed VF 

In this section, we compare theoretical VF in §[| with the observed one at z = 0©. We consider 
four specific cosmological models listed in Table |. The first three models (SC, OC and LC) use the 
scale invariant [n = 1) power spectrum of the cold dark matter (CDM) model in RefJ-3, and the last 
model (LP) uses that of the primeval isocurvature baryonciP (PIB) model in Ref. O) with the slope 
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Table I. Model parameters 

Model P(k) 7T \o h (T 1 n B fi 2 

SC CDM 1 05 L3 0.016 

OC CDM 0.3 0.7 0.52 0.016 

LC CDM 0.2 0.8 0.8 1.0 0.016 

LP PIB 0.2 0.8 0.8 1.0 0.128 



n = —1.15 at the large wavenumber. The Hubble constant h in each model is chosen so that the age 
of the universe be longer thari_L3 GyrEi) for given {Qq, Xq). The bias parameters b adopts the COBE 
DMR 2yr data normalization and the baryon density parameter fi-Q in the CDM models is from 
the 



as 



78) 



big-bang nucleosynthesis constraint W . The r parameter in the CDM model is empirically fitted 

(4-1) 



where Tq 



r = Q h{T /2.7K)- z exp[-i? B (l + V2h/f2 )] 
2.726K is the temperature of the cosmic microwave background. 



4> obs 

(Shimasaku 1993) 




v /(km/s) 



Fig. 1. PSVF rffiq.[3-6|) and KSVF (Eq.[3-7]) in the cosmological models listed in Table |. Crosses show the observed 
VF in Ref.E2l>. 

In Fig.l, we plot the theoretical VF's (PSVF and KSVF) for these models at z = 0. In Q + A = 1 
case, we used the fitting formulae for 5 c o_and $ v o in Eqs.( CT9 ) and ( |C-33| ). Plotted in crosses are 
the observed VF compiled by Shimasaku EHP. He multiplied the observed central velocity dispersion of 
stars in E and SO galaxies by a 1.5 1 ' 2 factor to take account of the difference between the distributions 
of luminous matter and isothermal dark halos. Since the error bars in the raw data represent the 
statistical ones only, we lengthened the error bars of galaxy VF by the 1.5 1 / 2 factor into the smaller 
velocity direction. However, it should be cautioned that this is an over-modification, especially at the 
low velocity end, because S galaxies do not have the factor. The cluster VF is translated from the 
observed temperature function of Ref.E3^ through the isothermal /3-model assuming (3 = 1. 

At least the VF observation should give a lower bound of the theoretical VF because the observed 
VF counts only luminous objects. In Fig.l, model SC fits the galaxy VF fairly well, but predicts too 
many large scale objects. Model OC fails on all scales due to the small value of the normalization ft -1 
from the COBE data. Models LP and LC are very good on cluster scale but significantly underpredict 
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galactic size objects, except for the KSVF in model LP. Shimasakuil) also compared PSVF with 
the observation, and concluded that no COBE normalized models are consistent with VF observation 
within PSVF. However, as noted in §3.3, PSVF has a flaw: PSVF neglects the accumulation of objects 
which formed earlier. Although KSVF and PSVF are not so different within the CDM models, they 
differ significantly in model LP. This originates from the fact that the PIB models have more powers 
on small scales than CDM models have: since smaller objects survive longer, the accumulation in PIB 
models is more remarkable than that in CDM models, especially on small scales. 

Nearly on all scales, KSVF in model LP seems consistent with observation if we allow the error 
bars to exclude the 1.5 1 / 2 factor. However, it should be borne in mind that there underlies a lot of 
uncertainties both in theory and observation of VF. On the observational side: (1) The 1.5 1//2 factor is 
quite generous and may represent the error considerably larger than the actual one. (2) Although the 
galaxies in the Virgo cluster is excluded in the ShimasakutJ data, it is possible that the galaxies in 
small groups are counted as individual objects in that data. From the construction of theoretical VF's, 
a group of galaxies should be counted as one object and the galaxies inside should be neglected. 13) 
The assumption of (3 = 1 for cluster VF may not be justified because of the /3-discrepancy problem EiP . 
(4) The observed VF may not represent the cosmic mean. On the theoretical side: (5) KSVF has some 
problems in counting argument 113. (6) We have not taken any dissipative processes into account in 
the dynamics of virialization, both in PSVF and KSVF. That is, we have related mass with velocity 
(Eq. |3-5|| ) only from a virial analysis and neglected a possible dissipation of energy from the collapsing 
object. 

Here let us discuss the last point (6) at some length. If one takes these dissipative processes into 
account, the velocity may increase because of the contraction to smaller radii. Based on the arguments 
of cooling and dynamical time scales H*' El 1 , these effects are expected to be unimportant on cluster 
scales. On galactic scales where the cooling is likely to be effective, theoretical VF will shift to the 
right and should become steeper as compared with the dissipationless case. This line of arguments 
might make the model LC consistent with the observed VF. We do not estimate here quantitatively to 
what extent the dissipation modifies the theoretical VF, but only note that the degree of dissipative 
effects in the galaxy formation can be measured by comparing the observed VF with the theoretically 
best VF without dissipation. 

In considering GL statistics, it is desirable to use VF's which are theoretically well-motivated 
and consistent with the observed VF. In particular, the evolution of the VF's is important in the 
calculation below. If one adopts the observed VF, however, one has to assume an ad-hoc number 
evolution. Given those uncertainties and problems, we have decided to apply theoretical and observed 
VF's in computing the GL probabilities on an equal basis. Further comparison between the theory 
and observation of VF should wait for future VF observation at higher redshifts. 

In Fig. 2 we plot the number evolution of PSVF and KSVF for the cosmological models in Table 
|. The decrease of small scale objects is slower in KSVF than that in PSVF, especially in models LC 
and LP. This is because the accumulation of small objects becomes less significant for higher redshifts 
in KSVF. The two VF's become almost identical at the redshift of ~ 4. Recently, Mo & FukugitatS 
compared PSVF with the abundance of giant galaxies at z = 3, and conclude that the spatially flat 
CDM models with Ao ~ 0.7 are favorable. However, their conclusion seems premature in that the 
model is already inconsistent with observation at z = (Fig.l). 



4.2. Comparison with the HST snapshot survey 



elL case, and compare the 



|),g^),|87),|3o|)_ Various uncertainties 



Let us now calculate the GL probability in the £ c = (SIS 
result with the observed GL frequency in the HST snapshot surveyl 
including the non-zero ii)re_raiiius are examined in the next section. Although there exist different 
surveys in optical bandit' E^'lHiJ, we do not mix these data to obtain the larger number of samples^ 
because we are mainly concerned with how the results are affected by various uncertainties described 
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Fig. 2. Evolution of KSVF (solid) and PSVF (dotted) for the cosmological models listed in Table |. SchVF without 
1.5 1 / 2 is also shown in dot-dashed curve for reference. From upper to lower curves in high velocity end, z = 0, 0.5, 
1, 1.5, 2, 2.5. 



in §4.1 and §g. 

Using the redshifts and V-band magnitudes of all the 502 QSO's in the HST survey, we calculate 
the expected number of lensed QSO's in the survey for various cosmological parameters as 



502 

nT P (8):=Y,P*(0,S u r n , Zl ) (4-2) 
i=l 



with Eqs. (EQJ ) and ( ^2l|) , where S t and Zi are the observed V-band flux and redshift of the i-th QSO 



in the survey^ The HST survey can detect multiple images if 9 > 0.1" and if the brightness ratio is 
smaller thanEP 

r„ = min[19(#/arcsec) - 85 , Si/S lim , 40] , (4-3) 
where Su m is the faintest flux that the HST can detect (the corresponding apparent magnitude is 



miim = 22 in V-band). We assume B — V = 0.2I53 when the V-band magnitudes are inserted into the 
B-band QSO luminosity function Eq.( [2-22| ). Theluminosities Li of QSO's are calculated from Si by 
Eq.(2-19) with (J?0)^o) = (1)0) and 7 = 0.5l^P'BlF. We use the standard angular diameter distance 
(q = 1; filled beam) for Da(z, z s ) in Eq.(2T5). 

Fig. 3a shows the prediction based on SchVF, while in Fig. 3b we plot the PSVF and KSVF 
predictions for the model parameters listed in Table||. In models OC and LC, the KSVF and PSVF 
predictions are almost indistinguishable (see Fig.l). Shown in the histogram is the distribution of the 
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Fig. 3. Image separation distribution of the expected number of lensed QSO's in the HST survey, (a): predictions from 
SchVF. (b): predictions from PSVF and KSVF. Model parameters in (b) are listed in Table |. Histogram shows the 
observed distribution. SIS lens model (£ c = 0) and the filled-beam distance (5 = 1) are used. 



observed 6 GL candidates in the HST survey. These are 1208+1011, 1413+117, 1115+080, 0142-100, 
0957+561, 1120+0154, and the image separation angles are 0.47", 1.22", 2.0", 2.2", 5.7" and 6.6", 
respectively. Despite the high resolution imaging of the HST Planetary Camera, 1208+1011 is the 
only sub-arcsecond GL candidate found by the survey©'©. 

Note that the SchVF can never predict the large separation lensing events such as 0957+561 and 
1120+0154, even if the 1.5 1 / 2 factor is included. Usually these large separation-events are neglected as 
statistical flukes in the literatures of GL statistics based on SchVF. Maoz & RixEP, for example, neglects 
0957+561 because it is not the lensing by a single galaxy but the lensing galaxy may be embedded in 
a clusterO'. They also neglect 1120+0154 because the images might be a physical pair (however, see 
Refs.t^'E^). However, we believe that it is important to consider these large separation lensing events 
when the PSVF or KSVF predictions are compared with observation, because the theoretical VF's in 
principle count all the gravitationally bound dark objects as well as luminous objects. The apparent 
absence of any obvious lensing object in 1120+0154 may be due to such non-luminous massive objects. 

It is also interesting to note that the predictions from observed VF and theoretical VF have 
opposite dependence on Qq. This is because theoretical VF naturally decreases with J?q (low density 
universe has fewer lenses). On the contrary, since SchVF does not at all depend on the cosmological 
parameters but is fixed by observation, the Qq and Ao dependence of SchVF prediction enters only 
through the cosmological geometry (distance formulae and the Hubble constant). In the case of PSVF 
and KSVF predictions, the Ao dependence comes mainly from the COBE normalized bias parameter 
b. 

Fig.3a shows that the SchVF prediction apparently prefers A = models to highly yl-dominated 
modelsB '!)'&■!). In Fig. 3b, model SC predicts too many lensing events in all separation range®. 
Again model OC is not viable: the COBE normalized density fluctuation amplitude in this model 
is too small to accommodate sufficient lensing objects. The most likely values of Qq in open universes 
should lie between 0.3 and 1. Model LC seems relatively good, while model LP is good at large 
separations but too much overpredicts the small separation events. We tabulate in Table II the 
expected number n exp := / n^ p {9)d9 of lensed sources for these models, and their \ 2 significance 
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level 0) 

r (n bin /2, X 2 /2) 2 _^« bs -n cxp ^ 2 
r(n bin /2) ' X ' ^ * 



Q== , X 2 := E V l exp' ^ , (44) 



where n° bs and n^ xp are the observed and expected number of lensed sources in the i-th bin. We set 
the bin width to be Z$bin = 1" as shown in Fig. 3, and calculate n exp and Q within the separation 
range 0" < 8 < 4" for SchVF predictions, and within 0" < 8 < 10" for PSVF and KSVF predictions 
(see the first row of the table denoted as "standard"). 

We note that, if the 1.5 1 ' 2 factor is excluded, the A = models in the SchVF prediction have 
very low significance level, because these models do not predict sufficient number of lensed QSO's 
around 8 ~ 2" (1115+080 and 0142-100). If the 1.5 1 / 2 factor is included, the values of Q in A = 
models are larger than those in A 7^ models only relatively: no models in the SchVF prediction can 
fit the observed distribution satisfactorily. The maximum likelihood analysis zB is unable to show this 
fact, where only the relative values of the likelihood among the modelsare calculated. This is why 
we do not use the likelihood analysis for the cosmological parameters o We caution, however, that 
the probabilities in the table are not statistically so meaningful because of the small number of lensed 
candidates. More reliable tests will be possible after the next-generation GL samples become available 
(see §|). 

As seen from Fig. 3 and Table II, GL statistics based on theoretical VF's prefers spatially flat low 
density models with the non- vanishing cosmological constant Jjf COBE normalized), as opposed to the 
SchVF prediction (this is consistent with the results in Ref.E$ which focused on the large separation 
lensing statistics using the observed cluster mass function in Ref.E^P). In particular, the model LC 
matches the observation fairly well. However, this is originated from the fact that the VF itself in 
the model LC underpredict the number of galaxy scale lenses today (Fig.l). Although the KSVF in 
model LP is the only theoretical VF (among the four models) consistent with the VF observation in 
Fig.l, Fig. 3b shows that the KSVF prediction of the model LP is unfavorably larger than the observed 
events at small separations. 

As discussed in §4.1 the observed VF on galaxy scales may be affected by dissipative effects. 
However, it is not clear which the observed VF or the theoretical VF without dissipation is relevant 
for GL statistics: if the lensing objects cannot be modeled as isothermal, then it might well be that 
the observed velocity dispersion in their central parts is larger than the effective value of v inside the 
Einstein radius of the isothermal model. If this is the case, the observed VF will overestimate the 
deflection angle and so the GL cross section. In passing we note that, if the observed VF is affected 
by dissipative effects, and if theoretical VF's without dissipation is relevant for GL statistics, then the 
model LC may become consistent with both of the VF and GL observations. 



§5. Various uncertainties in GL statistics 



Within the SIS lens model and the standard distance formula (a = 1), we could not find a 
cosmological model which satisfies both GL and VF observations. In this section we consider various 
uncertainties (free parameters) in the predictions of GL statistics. The results are summarized in 
Table II 

5.1. Distance formula 

We have no knowledge at present as to what value of the smoothness parameter a in the distance 
formula (Eq. |2-10fl ) describes our universe best. First of all, it is important to realize that the distances 
appearing in the lens equation (Eq. f2-2fl ) and the GL probability formula (Eq. p-12]j ) has slightly 

*^ We do use the Kolmogorov-Smirnov test here because we want to compare the amplitude of the observed event 
rate with that of the theoretical predictions. This is also why we use the \ 2 t es t although the binning is very crude. 



13 



Table II. The expected number n cxp of lensed QSO's in the HST snapshot survey and the x 2 significance level Q 
(Eq.Q) of the models (see Table |). (a,b): SchVF predictions within the separation range 0" < 9 < 4"; (c): PSVF 
predictions within 0" < 9 < 10"; (d): KSVF predictions within 0" < 9 < 10"; are compared with the observation. 



(a) VF 


SchVF without 1.5 1/2 


(J?o, Ao) 


(1,0) 


(0.3,0) 


(0.2,0.8) 


0" < 9 < 4" (n obs = 4) 


^exp Q 




^cxp Q 


Standard (§4.2) 


2.5 < 0.1% 


3.7 < 0.1% 


10 2.1% 


Distance formula 51 (§5.1) 


1.4 < 0.1% 


2.9 < 0.1% 


7.2 0.5% 


Core radius b (§5.2) 


2.3 < 0.1% 


3.4 < 0.1% 


9.4 1.2% 


Magnification bias c (§5.3) 


1.7 < 0.1% 


1.9 < 0.1% 


5.3 0.5% 


Dust obscuration d (§5.4) 


1.8 < 0.1% 


2.7 < 0.1% 


7.8 4.2% 



(b) VF 


SchVF with 1.5 1/2 


( Qq, Ao) 


(1,0) 


(0.3,0) 


(0.2,0.8) 


0" < 9 < 4" (n obs = 4) 


^cxp Q 


^cxp Q 


^cxp Q 


Standard (§4.2) 


5.1 38% 


7.5 30% 


20 0.7% 


Distance formula 3 (§5.1) 


2.8 3.7% 


5.9 25% 


15 5.2% 


Core radius b (§5.2) 


4.4 26% 


6.5 26% 


18 1.5% 


Magnification bias c (§5.3) 


2.6 9.1% 


3.8 19% 


11 19% 


Dust obscuration d (§5.4) 


4.1 45% 


6.0 43% 


17 2.7% 



(c) VF 


PSVF 


Model 


SC 


OC 


LC 


LP 


0" < 9 < 10" (n obs = 6) 


^cxp Q 




^cxp Q 


^cxp Q 


Standard (§4.2) 


29 3.2% 


1.7 < 0.1% 


6.1 79% 


16 45% 


Distance formula 3 (§5.1) 


17 55% 


1.3 < 0.1% 


4.6 48% 


9.3 81% 


Core radius b (§5.2) 


18 44% 


1.6 < 0.1% 


5.4 55% 


11 73% 


Magnification bias c (§5.3) 


23 17% 


1.4 < 0.1% 


5.1 60% 


11 79% 


Dust obscuration' 1 (§5.4) 


27 6.5% 


1.3 < 0.1% 


5.7 79% 


12 79% 



(d) VF 


KSVF 


Model 


SC 


OC 


LC 


LP 


0" < 9 < 10" (n obs = 6) 


^cxp Q 


^exp Q 


^exp Q 


^cxp Q 


Standard (§4.2) 


32 1.3% 


1.6 < 0.1% 


6.0 77% 


21 18% 


Distance formula 3 (§5.1) 


18 43% 


1.2 < 0.1% 


4.5 47% 


14 58% 


Core radius b (§5.2) 


21 27% 


0.7 < 0.1% 


5.4 53% 


19 20% 


Magnification bias c (§5.3) 


22 23% 


1.0 < 0.1% 


4.2 40% 


13 62% 


Dust obscuration' 1 (§5.4) 


29 3.4% 


1.2 < 0.1% 


5.6 77% 


19 33% 



a Self-consistent formula with a — is used. 

b Cc = 0.2(u/200kms- 1 ) 3 /i- 1 kpc. 

c Magnification of only the brighter image is used. 

d TcV = 1, £d = O.S/i^kpc. 



different meanings. First let us discuss why the filled beam (a = 1) distances are used in EaMl%). 

Actually Eq.(2-12) has two different derivations. One is the "self-consistent" formulaLiP which 
uses a = lin the GL probability formula (as in Eq.[ 2-12| ), and the other is the "optical-depth" 
formulae^' c3> which do not distinguish the distances in Eqs.( |2-2j ) and ( [2-12 ), and use the same value 
of a in both the equations. In the self-consistent formula, the lensing cross section is defined in the 
source plane (as in Eq.p-7|j) and the GL probability is defined by the ratio of the total lensing cross 



section to the whole area of survey regions in the source plane. In this case, the distances in Eq.(2-12) 
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measure the proper area in the lens and source planes subtended by the solid angle of the regions. Thus, 
if one observes GL frequency in the whole Am steradians in the sky, the filled beam distances should be 
used in Eq.( |2T2] ) because (we assume) the universe is on average homogeneous. On the other hand, in 
the optical-depth formula, the lensing cross section is defined in the lens plane and the GL probability 
is defined by the mean number of deflections along line-of-sight. Since D\{z s ) / Da{z s ) < 1, the optical 
depth formula is always smaller than the self-consistent formula©' |) and the difference becomes larger 
for smaller 5, larger z B , £2$ and Ao- For example, when 5 = 0, (i?,A) = (1,0), and z s = 2, the two 
GL probabilities differ by a factor of 0.645. We adopt the self-consistent formula because, as pointed 
out in Ref. Elf, the line-of-sight cannot be treated as a random variable but the source position on 
the source plane can be. Although the GL observation data used in §4.2 do not coye r th e whole 4tt 
steradians, the survey region is wide enough to assume the statistical homogeneity!^ '1129). 

On the other hand, the distances appearing in the lens equation (Eq, |2-2fl ) measure the p rope r 
lengths in the lens and source planes on scales of the typical image separation. It is shown (SasakiE21P) 



that the choice of 5 in Eq.(2-2) depends on how clumpy one imagines the universe is on th at s cale. 
So 5 in Eq. (|2-2| ) is a statistical quantity and may depend on the directions in the sky. SasakiliM* also 
argues that, when one calculates the deflection angle a = (d/d^)ip, the density profile of lensing objects 



should be subtracted by ap for consistency because the lens potential (Eq.[2-l|) is a perturbative 



quantity against the homogeneous background. In the case of strong lensing, however, the subtraction 
of ap makes little difference: if we assume that the lensing object has virialized before the light 
passes by it, the density of the lens is at least i? v times larger than the background when the light 
is deflected. Since the deflection angle is roughly proportional to the mass of the lens, and since 
i? v > 187T 2 (Eq. |]G18 |), the fractional change in the deflection angle due to ap is at most a factor of 



10 , which is observationally negligible. 

To be pedantic, setting a = 1 in Eq. Q2-2j ) is inconsistent because the strong GL cannot occur in 
the 5 = 1, exactly homogeneous universe. As 5 decreases from unity, the GL probability (without 
magnification bias) increases through the geometrical effect: because of the decreasing focusing, the 
proper area of the GL cross section increases. The effect of 5 is larger for larger Qq or Ao universes. 
This is because 5 in the Ricci focusing term (Eq. [p-10| ) is multiplied by Q(z) and larger Ao makes 
Q(z) approach unity faster at higher redshifts. On the contrary, if the magnification bias is taken into 
account, the GL probability decreases as 5 decreases, because of the factor p, := [Da(z s ) / Di(z s )] 2 in 
Eqs.( |2-21 ) and ( p-20| ). However, this decrease of GL probability with 5 is rather artificial: intuitively, 



as 5 decreases, the universe becomes clumpy so the number of lensing objects and the GL probability 
should increase. It may be that the number density of lensing objects should be multiplied by (1 — 5) 
to take this effect into account. 

Fig. 4 is the same as Fig. 3b except that 5 = (extreme case). Fig. 4a shows the PSVF and KSVF 
predictions from the self-consistent formula, while Fig. 4b uses the optical depth formula. Because of 
the decrease of the magnification bias mentioned above, the GL predictions decrease as compared with 
Fig. 3. In our opinion, previous work (e.g., Ref.EP) might overestimate the magnification bias because 
of the absence of the factor p in Eq.( [2-2i ): what is essential to the magnification bias is the relative 
magnification p(y)/p of a lensed source to unlensed sources. The x 2 significance levels of the models 
are tabulated in the second row of Table II. In extremely inhomogeneous (5 = 0) universes the models 
SC and LP are viable. 

5.2. Core radius 

It is known© that a small but finite core radius can change the lensing probability by an order of 
magnitude: since the strength of a lens is determined by its central condensation of the surface density, 
the removal of the central singularity reduces the lensing cross section significantly [oc (1 — x^ 3 ) 3 ; 
Eq.( |2-8] )1. However, as noted in Ref. eP, the inclusion of core radius increases the magnification bias: 
the significant decrease in the GL probability with the core radius is compensated to some extent 
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Fig. 4. Same as Fig. 3b, except that the smoothness parameter (Eq.[M0|) is a = 0. (a) self-consistent formula (Eq.[2-15]) 
with the empty-beam distance (a — 0) is used; (b): optical-depth formula with the empty-beam distance (d = 0) is 
used in stead of Eq.(2-15). 



by the increase of the magnification bias. The probability even increases with core radius when the 
source luminosity is greater than ~ 20L*. 

From Eq.( |2-3D , the scaled core radius x c is written as 



0.72 



100 km s" 



Dos 
Dls 



hD 



OL 



Gpc 



-i 



On galaxy scales the core radius is rather small. Fukugita & TurnerEf estimate £ c = 324 h pc at 



kpc 
i 



(5-1) 



v = 225 kms from the observation of E galaxies in Ref \ W2 v whichcor responds to x c ~ 0.12 for 
a lens redshift of z ~ 0.5. From the same data, Krauss & WhiteliS 3 ]) obtains £ c = 557/i _1 pc at 
v = 306 kms -1 (x c ~ 0.14) and £ c oc L ' 73 , though the scatter in the data is very large. Using the 
Faber-Jackson relation, one obtains p ~ 2.9 (see Eq.| 2-16| ). Wallington & NarayanElP argues that £ c 



must be smaller than about 200 h pc at v = 225 km s (x c ~ 0.08) so that the observed even number 
of lensed images is explained by the demagnification of the central image. On cluster seal 



pore 



radius is typically £ c ~ 20-30/i -1 kpc at v = 1000 kms -1 which corresponds to x c ~ 0.4-0.6^'^3)'li23). 
These data suggest that p ~ 3 and £ c ~ 0.2/i _1 kpc at v = 200kms _1 . 

Fig.5 plots the same quantities as in Fig. 3b except that the finite core radius is included and 
only the KSVF predictions are shown. We examined p = 1,2 and 3 (Eq. pT6fl ), in Figs. 5a, b, and 
c, respectively. From upper to lower of the five curves for each model, the core radius is ^ c = 
0, 0.1, 0.2, 0.3, 0.4/i _1 kpc at v = 200 kms -1 . When p < 2, the GL probability decreases mainly at 
small 6 and asymptotes to the SIS (£ c = 0) case at large 9, simply because x c oc v p ~ 2 (Eq. [p-6| ). When 
p > 2 the converse behaviour occurs, and when p = 2, it decreases uniformly. We note that the effect 
of the core radius is slightly smaller in larger Ao universes. This is because the scale length (Eq. |2-3j ]) 
increases with increasing Ao due to the cosmological geometry: larger means that the light traverses 
farther away from the lensing object, so the core radius looks effectively smaller. For the core radius 
of £ c = 0.2(?;/200kms -1 ) 3 /i -1 kpc, the x 2 significance levels of the models are summarized in the third 
row of Table II. 

The above estimates for the core radius of early-type galaxies are based on pre-HST data. But 
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Fig. 5. Same as Fig. 3b, except that the finite core radius £ c is included and only the KSVF predictions are 
shown, (a): £ c = f c »(v/200 kms" 1 ); (b): £ c = £ c »(v/200 kms -1 ) 2 ; (c): £ c = ^(v/200 kms" 1 ) 3 ; with f c , = 
0, 0.1, 0.2, 0.3, 0Ah~ kpc from upper to lower curves. 



recent observation of these galaxies by HST 1123 >[L2§) shows that they have nearly singular cores. More 



over, Ref.L^ZP argues that the inclusion of finite core radius should increase the velocity dispersion of 
dark matter for consistency. Based on these arguments, the effect of the core radius may be smaller 
than that given in Table II in practice. 

5.3. Magnification bias 

The magnification bias (Eqs. p-21 j and | 2-2C( |) contains which is defined in §2.1 as the total 



magnification of all the images. However, depending on the properties of a QSO sample or the GL 
configuration, fi(y) in the bias factor should be interpreted as the magnification of only the brighter 
image or the fainter image among the outer two images. If individual QSO's in a sample are not 
examined so closely whether they are lensed or not, the magnification of only the fainter image should 
be usedt^ because the fainter image should be bright enough to be recognized as one of the multiple 
images. On the contrary, if one examines the QSO's closely enough to search for the second image 
by a follow-up observation (as in the HST snapshot survey), the magnification of the brighter image 
or the total magnification should be usedl 108 P, Which magnification - brighter image or the total - 
should be used may depend on the image separation angle: if one searches for lensed sources of small 
separation angles, then the total magnification may be relevant, because it is likely that the brightness 
of a lensed source with a small separation is recognized as the total brightness of all the images. 

Fig. 6 shows how different choices for the magnification bias affect the GL prediction for the HST 
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Fig. 6. Same as Fig. 3, except that fj,(y) in the magnification bias (2-21) is defined as the magnification of the brighter 
image (the middle curves) or the fainter image (lower curves) of the outer two images. The upper curves use the 
total magnification of all the images for fj,(y) as in Fig3. 



survey. The upper curves for each model use the total magnification for n(y), while the middle 
and lower curves use the magnification of brighter and fainter image, respectively. Only the KSVF 
predictions and SchVF predictions without the 1.5 1 / 2 factor are shown. For the middle curves, the 
X 2 significance levels of the models are shown in the fourth row of Table II. Because the effect of the 
magnification bias is very large in the GL prediction, the different choices of the bias factor changes 
the result significantly. 

The uncertainties in the QSO luminosity function (Eq.[ 2-2^| ), a = 3.79 ± 0.15, = 3.15 ± 



0.10 and M * = -20.91 ± 0.25 + 51ogfc, change the GL prediction by ~ ±40%, ±40% and ±20%, 
respectively, while (3 = 1.44 ±0.20 does not cause significant changes because QSO's brighter than L* 
are prefe rent ially selected in the HST snapshot survey in order to achieve high efficiency of the lens 



detection 



109) 



5.4. Dust obscuration 

When the light from a lensed source passes by the lensing galaxy, the lensed image may be obscured 
by the dust ingredients in the galaxy. For example, the scale length (Eq. 

£* ~ 1.4/i _1 kpc 



100 km s" 



D LS hD OL 



is about 2/i _1 kpc for a lensing galaxy of v = 200 km s -1 and z = 0.5, which may be well inside the 
baryonic parts of the galaxy. The dust obscuration reduces the GL probability in two ways: i) due 
to the dust demagnification, the magnification bias decreases; ii) brightness ratio of the images is 
amplified and so it becomes hard to detect the second image, because the fainter image, which is 
nearer to the lens center, is more appreciably demagnified than thebrighter image. 

Let us estimate this effect quantitatively by a simple model E^'B'Ef. Since elliptical galaxies 
dominate over spiral galaxies in GL statistics (§3.1), the former is considered here as lensing galaxies. 
Let t(£) be the (2-dimensional) optical depth profile due to obscuration by dusts at the impact 
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parameter £ from the center of the galaxy. We assume the following form of t(£): 

r(0 = r c /[l + (£/£ d ) 2 ] (5-3) 

where r c is the optical depth at the center and £d is the dust core radius. The light ray which passes 
by the lensing galaxy at £ is demagnified by a factor e~ T ^\ Accordingly we have to redefine the 
total magnification fi(y) and the brightness ratio r(y) in Eq. fl2-21 ) by replacing fi p (x) in §2.1 with 
fj, p (x)e~ T ^* x \ andrecalculate the magnification bias. The optical depth is sensitive to the wavelength 
A. Following Ref. c£> , we use the fitting formula of Ref.fci5P: 

r c (A) ~r cV X(x)/3.2, (54) 

where x := (A/^tm) - ._7cV is the optical depth in V-band (xy = 1.8), and the functional form of 
X(x) is found in Ref. Ej). If QSO is observed in the wavelength A Q b s and the lensing galaxy is at 
red shift * then the demagnification factor should be evaluated at the wavelength A bs/(1 + z )- Wise 



& SilvaHiiP obtain 

T cV ~ 1 comparing their dust model with the data of 52 elliptical galaxies. We 
estimate £d ~ 0.3/i _1 kpc, assuming that £d is 0.1 times the effective radius of the de Vaucouleurs 
profilec2P. To be more realistic, r c y and £d should be dependent on the velocity dispersion v (the sizes 
of galaxies) and on z (the evolution of galaxies). However, we set them constant for simplicity because 
we aim at examining the extent to which the dust obscuration affects the GL probability. 



( d = 0.3h-'kpc (a) 



SchVF (CX3.0) 

••- (0.2,0.8) 

wilhoul 1.5 1/2 




2 4 6 8 

6 /arcsec 



Fig. 7. Same as Fig. 3, except that the dust obscuration effect by lensing galaxies is included, and only KSVF predictions 
and SchVF predictions without the 1 .5 1//2 fa ctor are shown. From upper to lower curves, the optical depth of the 



dust obscuration is r c v = 0, 1, 2, 3, 4 (see Eq.[5-3|) 



Fig. 7 shows the dust obscured GL predictions on the HST survey. Only KSVF predictions and 
SchVF predictions without the 1.5 1 / 2 factor are shown. Among the five curves for each model, the 
optical depth is r c y = 0, 1, 2, 3, 4 from the upper to lower curves, and £d is set to 0.3/i _1 kpc. The dust 
obscuration effect is considerably large at the small separation angles. We note, however, that the 
effect becomes small if one includes a finite core radius, because the brightness ratio becomes small. 
For the "realistic" value r c y = 1, the x 2 significance levels of the models are presented in the fifth row 
of Table II. 

Let us note general features of the effect of the dust obscuration. Larger Ao universes have smaller 
effect of the dust obscuration because the light rays are likely to pass by farther from the lens center, 
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as noted in §5.2. If £d increases with v less rapidly than oc v 2 , the GL probability decreases mainly at 
small separation and asymptotes to the no-dust curve at large separation, and the converse is true if 
£d increases more rapidly than v 2 , for exactly the same reason as in Fig. 5. If r c v increases with z due 
to the star formation activities at higher redshiftst^, lensed images with small separations are more 



in 



obscured than those with large separations, because t 
objects of high redshifts (see Eq. |2-6| |). Observation 
distributed as pdust(^) oc r~ 
obscured than in Fig. 7. 



re former are likely to be produced by lensing 
),|ii2p suggest that the dusts are more diffusely 
rather than as Eq.(5-3). In this case large separation images are more 



§6. Predictions for future surveys 



With the limited statistics of the current data, the constraints on the parameters are not yet 
decisive. We expect, however, that a much larger number of sources will be homogeneously sampled 
in one systematic lens survey in the near future. For example, the Sloan Digital Sky Survey (S PSS ) 
plans a spectroscopic survey of 10 5 QSO's over tt steradians brighter than 19 magnitude in g-bandlH^. 
In this section, we will make predictions on future GL surveys, bearing mainly the SDSS in mind. 

Let no(8,S)d0 be the expected number of lensed QSO's with image separation 6 ~ 6 + d9 and 
observed flux brighter than S, within the solid angle i? sa in the sky. We assume that the survey 
examines all the QSO's in the sky homogeneously. Then, using the cumulative QSO luminosity 
function in §2.4, ng(9,S) can be calculated as 

n e (e,S) = Q s J dz—r 2 (z)<P Q (L,z)P*(6,S,r*,z), (6-1) 
Jo dz 

where L is related to S by Eq.( |2-19| ), r(z) := (1 + z)D\{z) is the proper motion distance from us to 
z, and Eqs.( |2-18|) and ( 2-20| ) are substituted for Pf . Since we use Eq. (|2-22|) for ^q(L) which assumes 



(J?,A) = (1,0) and 7 = 0.5, these values of the parameters are also used in Eq.( (2-19| ) and r(z). We 
set z max = 5c3, and use the same detection limit r* = 19(#/arcsec) a85 as that of the HST survey. 




Fig. 8. Image separation 8 distribution of the expected number of all the lensed QSO's per unit steradian in the sky 
(Eq.[3-1]). (a): predictions from SchVF. (b): predictions from KSVF. From lower to upper of the three curves for 
each model, the limiting magnitude is m — 18, 19, 20 in B-band. 
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Table III. The expected number ol lensed QSO's per unit steradian brighter than 18, 19 and 20 B-mag. (a): SchVF 
predictions within the separation range 0" < 6 < 5"; (b): PSVF and KSVF predictions within 0" < 6 < 5"; (c): 
PSVF and KSVF predictions within 5" < 6 < 10"; for models in Table @. We adopt £ c = 0.2(u/200kms _1 ) 3 ^" 1 kpc, 
t c v = 1 and £d = 0.3/i -1 kpc as the models of the core radius and the dust obscuration (see §5.2 and §5.4). 
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Fig. 8 and Table III show the predictions on future surveys: the expected number of lensed QSO's 
per unit steradian (/? sa = 1)- We adopt £ c = 0.2(u/200kms _1 ) 3 /i _1 kpc, r c v = 1 and £d = 0.3ft. kpc 
as the models of the core radius and the dust obscuration, which we regard as the realistic values 
(see §5.2 and 5.4). The standard angular diameter distance (a = 1) is used. In Fig. 8, the KSVF 
predictions and the SchVF predictions without the 1.5 1 / 2 factor are shown. Among the three curves 
for each model, the apparent limiting magnitudes corresponding to S are m = 18, 19, 20 in B-band 
from lower to upper curves. Table Ilia shows the SchVF predictions within the separation range 
0" < 8 < 5", while Table Illb and IIIc show the PSVF and KSVF predictions within 0" < 9 < 5" and 
5" < 8 < 10". The KSVF prediction is strongly dependent on SIq and the bias parameter b. Therefore, 
we believe that Fig. 8b will provide a promising test for one of these parameters if the other is given 
independently. 

We propose that the cosmological tests should be carried out at a large separation range so that 
various uncertainties do not enter: the dust obscuration may be ineffective at 9 ^ 3" (Fig. 7); theoretical 
VF is free from uncertainties concerning the formation rate or the dissipative processes (§4.1) at 9 £ 7" 
(corresponding to v ~ 500 kms -1 ; Eqs.[2-6 and |5T]) which nearly discriminate the galaxy scale and 
cluster scale; etc. Alternatively, the lensing events with 9 ^ 7" should provide information about VF 
and the dissipative effects on the dynamics of galaxy formation. 

In closing we give some cautions on Fig. 8. It may be that the predictions in the figure are 
overestimated because of the assumption that the QSO luminosity function of Boyle et al.@* is valid 
up to z = 3. The dependences on the core radius and the dust obscuration are similar to Figs. 5 and 7. 
The uncertainties in the values of Mq and change the predictions by ~ ±20% and ±10%. Those 
of a and do not cause a significant change. 
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§7. Summary and conclusion 



it GL statistics disfavors cosmological models dominated by the 
)'EliQ), in contrast to the other cosmological tests suggesting low 



To date, it is well known tl 
positive cosmological constant A 

density universe with non-zero ylli3)'EJ). In the light of this, we reexamined the current constraints on 
A from GL statistics, mainly paying attention to the role of VF in the calculation of GL probability. 

Most previous work on GL statistics uses SchVF which counts only luminous galaxies as lensing 
objects. Moreover constant comoving number of lenses (no evolution) is often assumed. To take 
account of dark lensing objects and their number evolution, VF is constructed theoretically from the 
PS theory in §||EiP'E3'. Actually the PS theory in its original form is not a perfect one in deriving VF, 
and we attempted in §3.3 to construct a more realistic VF (KSVF) from the hierarchical clustering 
model in Ref.E^. Then we compared these theoretical VF's with the observation of VF by Ref.E^ for 
COBE normalized CDM and PIB models in §4.1. There we found that KSVF in the model LP - the 
spatially fiat PIB model with {[2q, Aq) = (0.2,0.8) (Table |) - matches the observation best (Fig.l). 
We also discussed that the model LC - the CDM model with (0.2, 0.8) - may become viable if the 
observed VF on galaxy scales is affected by the dissipative processes. 

At first sight, one might expect that the inclusion of dark lenses and their evolution into GL 
statistics would lead to tighter limits on Ao- -However, when we compared the predictions of the 
theoretical VF's with the HST snapshot surveyed in §4.2, we saw that spatially flat models with non- 
zero A (model LC) are viable among the COBE normalized models (Fig. 3). This is because the PSVF 
in low density universes underpredicts the number of galaxy scale lenses (Fig.l). Though the KSVF 
in the model LP is consistent with the VF observation, the model predicts higher lensing frequency 
at small separations than observed. We could find no COBE normalized cosmological model which 
satisfies both VF and GL observations, within the SIS, a = 1 and no dust model. 

In §[|, we examined various uncertainties in the GL predictions systematically: the choice of the a 
parameter in the distance formula (Eq. |2-10| ), core radius £ c , magnification bias, and the obscuration of 
images by dusts in the lensing galaxies. All these uncertainties reduce the GL prediction considerably 
from the standard calculation within a range of the realistic values of free parameters. The core radius 
and the dust obscuration are effective mainly at large and small separation range, respectively (Figs. 5 
and 7), while the uncertainties in a and the magnification bias reduce the GL prediction uniformly 
on all separation range. Because of these uncertainties, it is premature to put strong constraints on 
cosmological models, including A, from the rather small number of observed lensed sources. When a 
large number of QSO's are homogeneously sampled, a more reliable comparison will become possible. 
See Table II for a summary. 

In §|6], we made GL predictions relevant for the next generation survey like the SDSS. We believe 
that the predictions from theoretical VF will provide a promising test for Qq and Ao, combining with 
the COBE normalization. We proposed that the comparison between the predictions and the data 
should be made at large separation so that the various uncertainties do not enter: £ 7" may be ideal 
for cosmological tests because the theoretical VF is expected to be reliable on large scales. Lensing at 
the intermediate range (3" ^ 9 £ 7") may provide information about the degree of dissipative effects 
on the galaxy formation. At the same time, the theory of VF needs improvements on the formation 
rate (§3.3) as well as on the dissipative effects. 
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Appendix A 

Specific solutions of the Dyer-Roeder equation 



In this appendix, the specific solutions of the Dyer-Roeder equation (Eq. |2-1C| ]) are presented in 
the cases a = 1 (filled beam) and a = (empty beam). We use the units c = Hq = 1, and the 
notations a := (1 + z)^ 1 and 



w := Q 

with wo := f^Q 1 — 1 (see also Eq. [[Td|] ) 
i) a = 1 (filled beam) 
i-i) A = 



WqCL 
WqO? 



(A = 0) 

(n + x = i) 



(Al) 



D 1 (z 1 ,z 2 ) 
i-ii) Q + A = 1 



2a 2 



a{ /2 (l + «;i) 1/2 (l + 2^ 2 ) - a\ /2 {l + w 2 ) l ' 2 {\ + 2 Wl 



1/2, 



Dl(zi,Z2) 



2a 2 



1/2; 



1 1 7 



2' 6 6 



1/2; 



°i F o. «. «; -^i - < F o' «> «; -^2 



1 1 7 



2' 6 6 



(A-2) 



(A3) 



where .F is the hypergeometric function of type (2,1). 
ii) a = (empty beam) 
ii-i) A = 



D (zi,z 2 ) 



1 



, 2af (1 + ^)^(2^-3) 
-2a2 /2 (l + w 2 ) 1/2 (2u>2 - 3) + 3Kr 1/2 (??i - »&) 



(A4) 



where n := arccosh(l + 2iu) when I? < 1 and r/ := arccos(l + 2iu) when Q > 1. 
ii-ii) i? + A = 1 



^0(^1,^2) 



5/2; 



1 5 11 



5/2; 



fll F 1 2' 6' ¥ ; ~ Wl ' " ° 2 F I o' «> — ; ~ w * 



1 5 11 



2' 6 6 



(A5) 



We do not use Eqs.( [A3| ) and ( A-5|) because the evaluation of the hypergeometric function is much 
more time-consuming than the direct numerical integration of Eq. fl2-10|) . 

Appendix B 

Analytic expressions for v\ and v 2 



From Eqs.(|2-6|) and ( 2-16 ), v\ and v 2 in Eqs.( [2-13" ) and (2-15) are determined by the equation: 

v 4 - X 2 v 2p - Y 2 = , (B-l) 

where 



x . = 1 A)sgc* Y D os 



4-7T D lD ls 



8vr D LS 



(B-2) 
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The analytic expressions for v\ and V2 in the cases p = 0, 1, 2, 3 and 4 are: 

p = 0: Vl = {x 2 + Y 2 ) 1 / A 

p = 1 : Wl = [X 2 + (X 4 + 4y 2 ) 1 /2]i/2/^2 
p = 2: Wl = VY(1 - X 2 )- 1/4 



where 



When p < 2, v-i = oo. 



p = 3: 



■V2 



X 



X 



1 2 /2 

1 COS — 7T — Z 

3 3 \3 



n l/2 



1 j_ 2 7 

— — cos Z 
3 3 



nl/2 



p = 4 : « 12 = (^X)- 1 / 2 !! T (1 - 4X 2 y 2 ) 1 / 2 ] 1 /4 



Z := - arccos 



1-|(X 2 Y) 2 



(B-3) 
(B-4) 
(B-5) 

(B-6) 

(B-7) 
(B-8) 



Appendix C 

Spherical collapse in Q + A = 1 universe 

Consider a local spherical region of radius r and mass M in an Q + A = 1 universe with A > 0. 
From the expansion equation: 

d 2 r M A 



dt 2 



one obtains 



F , , f 2M A 



n H r . 

- 2 3 



2M yl ,^ _1/2 



(C-2) 



when dr/dt > (we neglect the decaying mode). The sphere "turns around" at a radius r ta defined 
by 



MA 

E = r 

rta 6 



ta 



if 



C := — — < - 
s 6M 2 



[our £ corresponds to n/2 of Ref.l 114 P], Then Eq.( p-2| ) is rewritten astli 



Ht 



A 



(i + C) + C^ 



1-1/2 



(C-3) 
(C-4) 

(C-5) 



where y := r/r ta . On the other hand, from Eq,( |l-2[ ) the age of the global universe in ,!? + A = 1 models 
is 

(C-6) 



Ht = -(l2w)" 1/2 arccosh(l + 2w) , 
3 



where w is defined in Eq.( AT ), The sphere collapses to r = in twice the turn-around time. In 
reality, however, it will reach a virialized state through the violent relaxation 0) and thereafter r 
stays constant r v determined by I 11 4 



4C^- 2 (l + C)y v + l = 0. 



(C-7) 
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where and hereafter the subscript v indicates the virialization time. Eq.(|C-7|) has the solution 



2 + 2C 
3C 



1/2 



1 



cos < - 7T — - arccos 
I o o 



1 / 3C \ 3/2 
C \2 + 2CJ 



(C- 



in the ranee < C < 0.5. Let p := 3Af/(4vrr 3 ) be the mean density inside the sphere. The over density 
against the background is written as 

p y 6 Q 

Let us calculate the overdensity of an virialized object at the virialization time: 

Pv _ Wv_ 
Pv 2/vC ' 



(C-9) 



(C-10) 



In order to relate ( with w v , we equate the local time inside the sphere (Eq. 



with the global 

cosmic age (Eq. |C-6| ) at the virialization epoch. The former is often approximated by the collapse 
time: 



H v t v 



1/2 „l 



dx 



--(i + C) + C^ 



-1/2 



-\- 1/2 [AK{k) - Bn{u, k)\ 



(C-ll) 
(C-12) 



where K(k) and n(u, k) are the complete elliptic integrals of the first and third kind, respectively (see 
Ref. i3 for their definitions), and 



A := [1 + (1 - 2C) 1/2 ]C, B := 2(1 - 2() l/2 C, 
C := 12C" 1/2 [2 - C + 2(1 - 2C) 1/2 r 1/2 
^2 ._ 2-C-2(l-2C) 1 / 2 



k 2 



v :- 



2 -C + 2(1 -2C)V2 ' 

-r 2 [i-c-a-2c) 1/2 ] 2 . 



(C-13) 
(C-14) 

(C-15) 

(C-16) 



Note that Eq.(C-12) does not assume the spatial flatness J? + A = 1. From Eqs.( C-6 ) and (C-12), ( 
and w v are related as 

Wv = -{cosh[AK(k) - Bn{u, k)} - 1} . (C-17) 

Now we obtain the analytic formula for i? v in fi + A = 1 universe through the parameter £. Given a 
value of C, one can compute # v and Q v from Eqs.( |C-^ ), flCTOp , ( pT7 ) and ( A-l ). In practice, however, 
it is useful to write $ v as a function of Q w . For this purpose, we give the fitting formulae: 

V ~ 18vr 2 ^ a573 or (C-18) 

i? v ~ 18^ 2 (1 + 0.40929w°' 90524 ). (C-19) 

Eqs.( P-l^ ) and ( P-19| ) are accurate within 5% for Q v > 0.1 and 1.7% for i? v > 0.01, respectively. 

Next, let us calculate 5 C , the density contrast in the early universe extrapolated linearly to the 
virialization time. Since y -C 1 in the early stages, we expand the integrand of Eq.(C_5) in powers of 
x and keep the linear term only: 



/ r\V 2 ry _ l 
Ht= \j) Jo ^ 1+ 2 {1+C)x]dx 



a) ^ 3/2 [ 1 + ^( 1 + ^ 



(C-20) 
(C-21) 
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Similarly, w -C 1 in Eq.flC-6|) yields 



3 



Equating Eq.(C-21) with ( C-22 ), y is written iteratively as 



1/3 



1/3 



(C-22) 



(C-23) 



Substituting this into Eq.( |C-9[ ), the density contrast of the spherical region in the early universe is 



6=~(l + ()(- 



w 



1/3 



(C-24) 



which grows as oc w 1 ^ 3 oc n. in accord with the linear perturbation theorylliZP. The linear growth rate 



in Q + A = 1 universe is 



118) 



(C-25) 



1/3 

which is normalized as D/a — > w when a — > 0, with i* 1 being the hypergeometric function of type 
(2,1). Extrapolating Eq.(C-24) by Eq.( |C-25| ) until a v , one obtains 



1/3 



(C-26) 



One can compute <5 C and 17 V from Egs- flC-l"^ ), ( |C-26| ) and ( A-l ) through the parameter £. We give 
fitting formulae so that one can obtain 5 C directly from a given value of i? v : 



Sc ~ Ji(i2vr) 2 / 3 ^.00539 or 
20 



(C-27) 



5 C ~ — (12vr) 2 / 3 (l + 0.012299 log 10 Ar). 



(C-28) 



Eqs.( P-27p and ( P-28j ) are accurate within 2.4% for Q v > 0.1 and 0.1% for f2 v > 0.01, respectively. 
Of more practical use (in, e.g., the PS theory) are the following quantities 



■#vO := Pv/Po = #v/ctv 



(C-29) 



D . 3 „ /l . 11 \ ^ /wn\ 1/3 



^-^-FU^-f'-^JP+oCTJ (C ' 30) 

= ^ (12 ' r)2/3 i F G' 1 'T ; "'" , ') x( ''°- ) ' ( ° 31) 

where 

X(w v ) ~ f?-°- 215 or (C-32) 

IK) ~ 1 + (5.1066uv L0812 + 2.02 15U;- - 35396 )- 1 . (C-33) 

The accuracies of Eqs. (|C-32|) and (|C-33|) are within 1.8% for [2 V > 0.1 and 3.1% for [2 V > 0.01, 
respectively. The corresponding formulae in an open (A = 0) universe are found in Ref. ©. 
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